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Abstract 
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tH- ' We present a new technique developed for the measurement of the scintillation 

On . efficiency of nuclear recoils in solid scintillators. Using this technique we measured 

the quenching of the scintillation efficiency for the various recoiling nuclei in CaW04 
^^ ' crystals which are used in direct Dark Matter searches. 
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Introduction 



The direct detection and tlie understanding of the nature of Dark Mat- 
ter remain the outstanding challenges of present day physics and cosmology. 
Searches for Weakly Interacting Massive Particles (WIMPs) which may con- 
stitute the galactic Dark Matter are presently being carried out by a number 
of groups world-wide [1][2][3]. Direct detection experiments appear today as 
one of the most promising techniques to detect particle Dark Matter. All di- 
rect detection experiments rely on the basic idea that many WIMPs should 
pass through the Earth, making their detection possible. Interaction with or- 
dinary matter is expected via elastic scattering on nuclei which, for the spin- 
independent case, should be coherent with all nucleons in a nucleus resulting 
in a cross-section rising in proportion to A^ (where A is the mass number) 
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and therefore favoring heavy nuclei [4]. In addition a low energy threshold 
is an experimental requirement as the predicted energy spectrum falls expo- 
nentially with energy. Due to the extremely low WIMP -nucleus scattering 
cross section the ability to discriminate WIMP signals from the radioactive 
background is essential for this kind of experiments. 

Cryogenic detectors based on simultaneous detection of light and heat 
fulfill these requirements [5]. They utilize the fact that the light output from 
nuclear recoil events in the scintillator is reduced compared to the light ob- 
served from f3 and 7 - interactions ( "electron recoils" ) of the same energy. This 
is characterized by the Quenching Factor (QF) which is defined as the ratio 
of the light output from electron recoils and from nuclear recoils of the same 
energy. For an effective background suppression the quenching factors have 
to be known, especially because different nuclei present in the scintillator are 
expected to have different quenching factors [6] [7]. In contrast to the WIMPs, 
the residual neutron background comes mostly from elastic scattering off the 
lighter elements. Therefore the different quenching of different recoiling nuclei 
in the scintillator can be used for a discrimination of the neutron background. 

In order to measure the quenching factors of different recoiling nuclei in 
scintillators we have developed a new technique. It relies on the measurement 
of the light output produced when the scintillating crystal is bombarded by 
different ions. This method allows to measure quenching factors for a wide 
range of different elements and scintillators. 



1 Experimental technique 



Nuclear recoil quenching factors have been measured for a variety of scin- 
tillators such as Nal [8] [9], CaF2(Eu) [10], CsI(Tl) [11] with neutron scattering 
experiments. 

Instead of measuring the light produced by the recoiling nuclei of the tar- 
get material as in the neutron scattering experiments, the technique presented 
here measures the light emitted when single ions impinge onto the scintillating 
crystal. As this work was done in the framework of the CRESST ^ experi- 
ment, CaW04 was used as target material. Ions with fixed kinetic energy were 
produced with a Laser Desorption / Ionization ion source (LDI) in a time-of- 
flight mass spectrometer (TOF-MS). The CaW04 crystal was mounted at the 
end of the flight tube and the light was detected using a photomultiplier. An 
^^Fe source was used for the reference measurement. 
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1.1 TOF mass spectrometer with an LDI ion source 



Figure 1 shows the scheme of the time-of-flight mass spectrometer with 
the LDI ion source used for this measurement. The ion source was separated 
with a small gate valve from the main vacuum chamber of the spectrometer. 
This allowed to vent just the small compartment for the exchange of the target, 
with minimal interruption (^ 20 min) of operation. 

Six different target materials could be fixed on the target holder, using sil- 
ver epoxy (Fig. 2). The target holder could be rotated and moved transversely 
to be able to reach the whole target with the laser. 

A shot from the pulsed UV laser of 337 nm wavelength desorbed and 
ionized the atoms. The shot duration was in the order of nanoseconds. The 
laser power could be controlled with a dielectric attenuator and the beam was 
focused on the sample disc with a lens of 200 mm focal length. At the sample 
target the laser power density was typically 10^-10^ W/cm^. The maximum 
acceleration voltage of the LDI ion source was 18keV. The ion beam was 
focused on the scintillator using an Einzel lens and X-Y deflection plates 
(Fig.l). 

The MS was instrumented with an additional ion deflector which allowed 
the selection of ions with a certain mass. It consisted of a pair of plates the 
electric field of which was turned on when the laser shot was fired. It deflected 
all ions out of the beam path until the electric field was switched off. Therefore, 
all lighter ions preceding the ions of interest were deflected. The high voltage 
was then switched off to allow the passage of the desired ions and then switched 
on again in order to block higher mass ions (see Fig. 3). 

The time-of-flight part used for the quenching factor measurements con- 
sisted of two ion reflectors (for details see [13]). Between these two reflectors a 
micro channel plate (MCP) could be inserted into the ion path at any moment 
of operation (Fig. 1). It was used to define the time window for switching off 
the voltage applied to the deflection plates to allow only the passage of the 
ion species of interest from a contaminated target material or from a multiele- 
ment target (see Fig. 3). The selected elements and target materials are listed 
in table 1. The time-of-flight measured with the MCP, which is roughly half of 
the one at the outlet of the mass spectrometer, was used to adjust the photon 
counting window appropriately for each element (see Sec. 1.2 and Fig. 5). 

The laser pulse was used as a start trigger and the signals of the MCP for 
the time-of-flight measurement were read out using a 500 MHz PCI transient 
digitizer (Fast Com Tec - TRP 250/500). 



1.2 PMT setup for photon counting 



The photon counting setup consisted of a small CaW04 crystal (5x5x5 mm"^) 
held in a teflon reflector and viewed by a photomultiplier tube that was placed 
in a housing which was mounted at a gate valve at the outlet of the mass 
spectrometer (Fig. 1 and Fig. 4). Ions reached the CaW04 crystal via a small 
hole (0 = 1 mm, length = 8 mm) in the reflector. In order to reduce the effect 
of trapped light in the crystal, the crystal surface facing the photomultiplier 
tube (PMT) was roughened ^ . Additionally, a thin layer of silicon grease was 
placed between the crystal and the PMT for optical coupling. The PMT was 
an Electron Tubes 9124B selected for a low rate of dark counts. A preamplifier 
followed by a 500 MHz digitizer was used for the read out of the PMT signal. 
The data acquisition was triggered by the laser. For each laser shot the number 
of photons was determined by counting the number of voltage pulses exceeding 
a given threshold amplitude in a fixed time window after the arrival of the 
ion type of interest (see Fig. 5). The number of registered photons for each 
shot was then accumulated in a histogram. Additionally, for each registered 
photon its arrival time was collected in a second histogram, which gave the 
light curve of the scintillator and the time-of-fiight of the ions (see Fig. 7). 

The reference measurement for the quenching factor was made with 5.9 keV 
X-rays from an ^^Fe source. The radioactive source could be installed in front 
of the collimator without any change of the setup. For the reference measure- 
ment the data taking was triggered by the PMT signal. The photon counting 
was performed as in the case of ions. 



2 Results and Discussion 



2. 1 MCP measurements 



The elements (H, O, Si, Ca, Cu, Y, Ag, Sm, W, Au) were chosen to 
smoothly cover the broad mass range of the elements in CaW04 (see table 1). 
For a proper placement of the time window for photon counting the arrival 
time of each ion type was measured with the MCP. As shown in Fig. 6 (left) the 
arrival times measured with the MCP follow the expected linear dependence 
on ^fA ^ . 
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2.2 PMT measurements 



The width of the counting time window was kept fixed (30, 40, 50 /xs) 
for all measured elements. Figure 7-left shows light curves of CaW04 obtained 
with beryllium and copper ions. The times of flight extracted from the onset of 
the light curve (Fig. 6 right) were roughly twice as long as the ones measured 
with the MCP, consistent with the longer flight distance. 

The sharp rise time (< 50 ns) of the recorded light curves conflrms the 
absence of electrical fields due to a charge up of the crystal or of its tefion 
holder during the measurement. Such charging up would slow down the ions 
before hitting the scintillator and therefore change their time of fiight and 
this in turn would broaden the rise of the light curve or shift the onset posi- 
tion of the light curve in consecutive measurements. A comparison of several 
subsequently measured light curves confirmed the absence of such charging 
up. 

Figure 7-right shows the biexponential fit of the measured light curve. 
The faster decay time determined this way is 1.08 /is while the longer one is 
7.96 /is. Measured decay times are in the agreement with recently reported 
values [12]. 



2.2.1 Reference measurement with X-rays 

The electron recoil reference measurement for the determination of the 
quenching factors was made with 5.9 keV X-rays from an ^^Fe source. Figure 8 
shows the spectrum measured with a 40 //s photon counting window. Simul- 
taneous measurements with three time windows (30, 40, 50 /is) were made in 
order to probe the dependence of the results on the width of the chosen time 
window. The spectra have been fitted using the fit procedure developed for 
ion spectra which will be discussed later. The results are presented in Table 2. 



2.2.2 Spectra measured with ions 

The laser power and the ion beam optics were adjusted such that the 
probability of multiple ion events (more than one ion impinging on the crystal 
surface per laser shot) was reduced as much as possible. The procedure is 
illustrated in figure 9. Laser power and focus were tuned that on average only 
each 15th laser shot gave an ion signal (Fig. 9 bottom). 

where A is the mass number of the ions, U is an electric potential difference, z 
charge state of the ions and D is the flight distance. 



2.2.3 Data analysis 

For each laser shot a certain number of photons was detected. The main 
contribution to the recorded spectra comes from laser shots not producing ion 
signals and thus measure the dark counts of the photomultiplier. The number 
of dark counts obeys a Poisson distribution with a small mean value. The 
number of the photons produced by single ion events ^ gives a spectrum with 
a Poisson distribution with a higher mean value. If multiple ion events exist 
additional Poissonians will appear with a mean value equal to a multiple of 
the mean of the first one. In the counting procedure the origin of the registered 
photons cannot be identified. Therefore, the background had to be accounted 
for in the fit function. 

For simplicity first only none (background) and one ion (signal) arriving per 
laser shot will be considered. In this case the probability to have k registered 
photons, Pt(k), is given by 

fc 
Pt{k) = P.on(0) ■ P{k,B) + P,on{l) E P{l,B)Ps{k-l,n ■ 1), (2) 

1=0 



where Pion{i) is the probability for i ions arriving per laser shot, P{k, B) is 
the probability of observing k background counts from a Poissonian distribu- 
tion with mean B and Ps{m,n) is the probability that m signal counts are 
observed when n is the average number of observed photons for single ion 
events. The term describing the ion signal sums the probabilities of all possi- 
ble contributions of background and signal counts. When up to A^ ions arrive 
per laser shot equation 2 generalizes to 

N k 

Pt{k) = P,,„(0) ■ P{k,B) + Y.P^on{i) E P{l.B)Ps{k-l,n ■ i) . (3) 

i=l 1=0 



Here, Pioni^), n and B are the parameters to be determine by the experiment. 
The background is described by a Poissonian with a scale factor Pjon(O). The 
mean value B of the background Poissonian is determined by a separate back- 
ground measurement for each time window. The ion signal can be fitted with 
a selectable number of multiple ion events to be taken into consideration. The 
scale factor Pioni'i) for each of the signal Poissonians i.e. the probability that a 
number i of ions arrives per laser shot and the average number n of detected 
photoelectrons per ion are the remaining free fit parameters. In figure 10 dif- 
ferent fits to the photon multiplicity spectra measured for H and W ions are 
shown. The dashed lines represent fitted spectra using the method described 
above. 



^ Due to the reduced laser power and applied beam defocusing only single and a 
smaller number of double ion events are occurring. 



As the peak positions of the single and double ion events of the hydrogen 
spectra are well separated from the background it was possible to check the 
validity of the background model used in the fit. Therefore, the three scale 
factors (Pion(O), Pjon(l), -Pion(2)), u and B were left as free fit parameters. 
Within the statistical errors the values obtained for B are in agreement with 
those from the separate background measurements (for the 30 /iS time widow 
determined B values are: 0.0180±0.8E-03 from the separate background mea- 
surement while 0.0176±0.12E-02 from the H spectrum.). To avoid the influence 
of the correlation between fit parameters n and B on the result of n, the B 
values determined by the separate background measurements were used in the 
analysis of all ion spectra. 

In the case of the hydrogen spectra (Fig. 10 left) there is a definite excess 
of measured events over the fitted line on the left hand side of the single ion 
peak ^ . Such a background could be either produced by a relatively small 
number of ions that reached the crystal surface with a reduced energy or by 
backscattered ions which deposit only part of their energy in the crystal. In 
order to address this question, the process of the ions impinging onto the crys- 
tal surface was simulated using the SRIM2003 simulation package (Stopping 
and Range of Ions in Matter version 2003 ) [14]. This package calculates the 
stopping power and the range of ions in matter using a full quantum mechani- 
cal treatment of ion-atom collisions. Here, only the effect of backscattered ions 
will be discussed. 

Figure 11 shows the simulated spectra of the energy deposition in the 
target (CaW04) crystal from the backscattered ions for several projectiles. 
The total number of simulated ions was the same for all elements (10000 
ions) and the numbers of backscattered ones are indicated in the figure. The 
energy spectra as well as the probabilities for the backscattering obtained from 
the simulation were implemented in the fit using the following procedure. An 
energy deposition with an energy lower than the nominal one will result in the 
production of a reduced number of photons obeying Poisson statistics within 
each energy bin. In order to consider energy losses due to the backscattering, 
the following fit function was used 

TV k ^ion 

Pt{k) = P^on{^)P{KB) + Y,P,on{^Y.P^^^^) [ Ps{k - l,ni^)PBs{Ed)dEd , (4) 



where PBs{Ed) is the probability density that due to the backscattering losses 
the energy E^ is deposited instead of the full energy of the ion Eion- PssiEd) 
has been calculated with SRIM2003. Fits including the effect of backscattering 



^ This effect was only observed in the case of the hydrogen spectra where the 
signal peak position is far from the background, whereas for the heavier ones the 
Poissonians overlap such that the effect can not be noticed. 



are shown in figure 10 as solid lines. Adding the backscattering effect gave a 
better agreement with the measured data in the case of the lighter elements. 
Therefore, the backscattering of the ions can be only a partial explanation of 
the excess events observed on the left side of the single ion peak. No change 
was observed for the heavier ones where the backscattering effect is only in 
the permil range. 

For the light ions like hydrogen another possibility for energy losses is 
conceivable. Ions hitting the collimator inner surface could be scattered with 
a shallow angle and therefore reach the crystal surface with a reduced energy. 
In order to address this issue, a test measurement was performed where a thin 
kapton tape was glued on top of the existing ion collimator hole (0=1 mm). A 
hole with a diameter of 0.7 mm was made in the center to prevent the ions from 
hitting the surface of the long collimator hole. This way, the number of the 
events that can scatter off the teflon surface should be considerably reduced. 
The resulting spectra showed a lower relative intensity of the excess events, 
but most probably due to a non perfect alignment of the used collimator, the 
effect was still partially observed. This indicates that the nature of the small 
excess background, observed in the case of H ions, in the region below the 
single ion peak can be at least partially explained, but further investigation is 
needed for a full understanding. Nevertheless, these excess events are only a 
very small fraction of the total counts and have a negligible influence on the 
number of photoelectrons per single ion hit resulting from the fit. 

In the case of heavier ions the fit was more sensitive to the number of 
Poisson distributions included in the fit due to the strong overlap of the dis- 
tributions (see Fig. 10). As previously discussed, via the consistent use of de- 
focusing of the beam for all ions the probability for multiple ion events was 
reduced to a very low level and therefore only single and a small contribution 
of double ion events had to be taken into account for the determination of n. 
The introduction of a triple ion contribution into a fit of the W data leads 
to a strong correlation of the fit parameters resulting in unphysical intensity 
relations (Pion(3) > Pion(2)) (see Fig. 10). Nevertheless, the resulting n agree 
within statistical errors. 

Technically it was not possible to obtain a beam with only and 1 ion 
arriving due to the sharp dependence on the laser beam power and defocusing. 
This effect is illustrated in figure 9, where even for the smallest laser power 
and the weak focusing a certain contribution of two ion signals is apparent. 
Therefore the single ion model is excluded. 

According to the previously discussed points, we decided to include the 
backscattering effect and to take into account only single and double ion con- 
tributions in the final fitting procedure. 



The ratio of the number of registered photons from the ion spectra, 
pe~{ion), and the reference measurement, pe~{^^Fe), normahzed (to IkeV), 
will give the quenching factor for that specific element: 

pe-(55Fe) ISkeV 

^^(^^") = ^Fi^^ ■ 5.9576 A:ey ' ^'^ 



where 5.9576 keV is the weighted energy of K-lines from ^^Mn [15]. 

The linearity of the light output of the used sample was studied in a 
separate measurement with various gamma sources. The result is shown in 
figure 12. The observed good linearity justifies the reference measurement at 
5.9 keV. 



2.2.4 Quenching factor results 

Figure 13 shows fitted spectra for oxygen, calcium and tungsten ions of 
18 keV measured with a 40yus time window. Mean numbers of photons for 
single ion events obtained from the fits to spectra of all elements, as well 
as numerical values calculated for the quenching factors are given in table 4. 
A graphical representation of the dependence of the quenching factors on 
the atomic mass of the element, for 40 /is time window, is given in figure 14. 
Within the 30 yus time window not all emitted light is collected (see Fig. 7) 
whereas practically no light is emitted more than 40 ^us after the impact. This 
explains the slightly different mean values of pe~ measured in the 30 ^s and 
40 //s time windows (see table 4). Practically identical values are obtained in 
the 40 yus and 50 yus time window measurements. Figure 14 also includes results 
from neutron scattering experiment at room temperature [16] and the values 
measured with a cryogenic detector at 7mK [3]. The presented data for oxy- 
gen and the neutron scattering result are in reasonable agreement. The low 
temperature values follow the systematic mass dependence of the room tem- 
perature results. This demonstrates the absence of a significant temperature 
dependence of the quenching factors in CaW04. 

The presented new method for the quenching factor measurements has 
several advantages compared to traditional neutron scattering experiments. 
The main advantage is that it is a simple "table-top" experiment. Measure- 
ments can be performed within a short period of time and a profuse choice 
of nuclei. A possible limitation may be due to surface effects causing a degra- 
dation of the light output for energy depositions close to the crystal surface. 
Therefore, the technique can be used in principle for any (solid) scintillation 
material, which does not suffer from a surface degradation of the scintillation 
efficiency. In the case of CaW04 crystals we could not observe any surface 
dependence. Strong arguments for this could be found in the comparison of 



the light yields from the ^^''Po alpha decays at the crystal surface and ^^°Po 
decays in the bulk of the crystal measured with cryogenic detectors (see [3, 
Fig. 7 ]). The light yield is the same for both groups of events, showing, at 
least for alpha particles, that the quenching factor is the same for interior and 
surface events. This supplements the argument in [5] that the absence of any 
splitting of the electron-photon band, where one has electrons mostly near 
the surface and photons mostly in the interior, demonstrates that bulk and 
surface events have the same quenching factor. 



3 Conclusion 



The quenching of the scintillation light from nuclear recoils is the key for 
an efficient background rejection in cold Dark Matter particle detectors that 
utilize scintillators. A new and powerful technique for the measurement of 
quenching factors was developed. It relies on the measurement of the number 
of photons produced when the scintillating crystal is bombarded by different 
ions. It uses a mass spectrometer to accelerate ions and measures, in single 
photon counting mode, the light which they produce when impinging onto a 
crystal. 

With the new technique presented here it was possible for the first time 
to determine the quenching factor of tungsten in CaW04. The measurement 
of quenching factors of different recoiling nuclei allows their identification and 
offers practically a multiple target choice within one absorber crystal. This 
provides additional sensitivity by discrimination of neutron background (O - 
recoils) and possible WIMP events (W- recoils) [3]. In the case of a positive 
Dark Matter signal the change in the recoil energy spectra of different absorber 
nuclei gives a unique signature and a powerful tool for further identification 
and verification of WIMPs properties. 
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Table 1 

List of selected elements with associated target material. 



Ion 


Target material 


Ion 


Target material 


H 


Stainless steel 


Y 


Yfoil 


Be 


Cu-Be foil 


Ag 


Ag foil 





CaO powder or CaW04 crystal 


Sm 


Sm foil 


Si 


Stainless steel 


W 


CuW foil and CaW04 crystal 


Ca 


CaO powder or CaW04 crystal 


Au 


Au foil 


Cu 


Cu-Be foil 







Table 2 

Fitted values of the numbers of photoelectrons for the ^^Fe peak. Three different 

time windows have been analyzed. The statistical errors {la) of the fits are quoted. 



Time windo^v 


Mean value 


30 ^s 


23.67 ±0.13 pe^ 


40;us 


24.08 ± 0.13 pe" 


50 /US 


24.22 ± 0.14 pe" 
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Table 3 

Comparison of fit parameters using different fitting procedures shown in figure 10. 
Tlie statistical errors (la) of the fits are quoted, x^ /dF is the chi squared per degrees 
of freedom and dF degrees of freedom. Discussion is given in the text. 



Elcm. 


Model 


xVdF 


dF 


n 


P«o„(l) 


P»o„{2) 


P.on(3) 


H 


Double ion events and 
backscattering effect 


57.2 


73 


32.96±0.11 


3968±64 


1117±37 




Triple ion events and 
backscattering effect 


39.0 


108 


33.00±0.11 


3963±64 


1149±37 


425±24 


Double ion events and 
NO backscattering effect 


169223 


71 


31.71±0.10 


3894±64 


1092±37 




W 


Double ion events and 
backscattering effect 


6.4 


3 


1.81±0.11 


1324±122 


611±126 




Triple ion events and 
backscattering effect 


2.0 


4 


1.79±0.14 


1551±122 


142±37 


361±80 


Double ion events and 
NO backscattering effect 


6.4 


3 


1.81±0.11 


1324±122 


612±126 
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Table 4 

Mean number of detected photoelectrons and calculated quenching factor values. 
The statistical errors (la) of the fits are quoted, x^ /dF is the chi squared per degrees 
of freedom, dF degrees of freedom and TW is the measurement time window. 



Element 


TW 


n 


P>on(l) 


Pzon(2) 


xVdF 


dF 


P»on{0) 


QF 


H 


30^s 


32.96 ± 0.11 


3968 it 64 


1117 it 36 


57.2 


73 


34079itl85 


2.17 ± 0.02 


40^s 


33.77 ± 0.11 


3965 it 64 


1110 it 37 


56.7 


75 


34062itl85 


2.15 ± 0.02 


50^s 


33.56 ± 0.11 


3933 it 64 


1105 it 37 


149.7 


75 


34135itl85 


2.18 ± 0.02 


O 


30^s 


5.16 it 0.07 


3082 ± 93 


2436 ± 86 


24.8 


12 


52011di230 


13.9 ± 0.2 
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Fig. 1. Scheme of the experimental setup. 
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Fig. 2. Target holder with six different materials. 
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Fig. 3. Application of deflection plates for the selection of desired ions. Top: TOF 
spectrum of a CuBe target without deflection measured with the MCP. Applying 
the proper selection parameters only ions of beryllium (middle) or copper (bottom) 
are reaching the detector (and also the PM setup at the end of the TOF path). 
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Fig. 4. Schematical view of the PM setup for photon counting (not to scale). The 
components are discussed in the text. 
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Fig. 5. Photoniultiplier signal from one laser shot. The data acquisition is triggered 
by the laser. After the time-of- flight the ion hits the crystal. In the figure voltage 
pulses from single photons created by the ion are shown. Ions are counted within a 
fixed time window which opens shortly before the ion arrives after the known time 
of flight. The width of the counting window was kept constant for all measurements. 




Fig. 6. Dependence of time-of-flight on \/A (where ^ is a mass number) measured 
with the MCP placed at half the distance to the PM tube (a) and with a CaW04 
crystal in the PM setup placed at the outlet of the mass spectrometer taking the 
onset of the light curve as the arrival time of the ion (b). The solid lines represent 
linear fits to the measurements with parameters as given in the boxes. 
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Fig. 7. Light curves for beryllium and copper (left) and bi-exponential fit of the light 
curve for beryllium (right). The fit parameters are indicated on the graph. The fast 
decay component is 1.08 ^s while the slow one is 7.96 /is. 
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Fig. 8. Photon multiplicity spectrum recorded with a time window of 40 /is while 
irradiating the CaW04 crystal with X-rays from ^^Fe source. The corresponding 
fitted Poissonian mean value is 24.08±0.13 pe" (4.04±0.02pe"/keV). 
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Fig. 9. Illustration of the effect of the beam defocusing on the shape of photon 
multiplicity spectra measured for hydrogen ions of 18keV. The total number of 
laser shots for the highest focusing is ~ 8500 (top) while for the other two spectra 
this number is doubled. Defocusing reduces the probability of multiple ion events 
and significantly prolongs the measurement time. Technically it was not possible to 
obtain exclusively single ion production. 
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Fig. 10. Comparison of different fits to photon multiplicity spectra. Two examples 
are shown: hydrogen, 18 keV (left) and tungsten, 18 keV (right). In the case of hydro- 
gen a higher laser intensity was chosen on purpose to produce a visible contribution 
of two and three ion events. For the tungsten case a lower laser intensity was chosen 
producing a small contribution of double ion and negligible contribution of triple ion 
events. The squares represent the measured data. The solid lines are likelihood fits 
including single and double ion events and backscattering of the ions. The dashed 
lines are fits without the backscattering effect while the dotted lines show the fits 
which include triple ion events and backscattering of the ions. The data point with 
zero photon counts is far above the range shown. The fit parameters are given in 
boxes and table 3. The discussion is given in the text. 
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Fig. 11. Reduced energy deposition due to backscattered ions. A total of 10000 
ions has been simulated with SRIM2003 package in each case. For each element 
the number of ion hits with a reduced energy deposition due to backscattering is 
indicated in the upper left corners of the graphs. 
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Fig. 12. Linearity of the light response of the CaW04 sample used for the measure- 
ments. The CaW04 crystal was irradiated with several 7 sources indicated on the 
graph and the photon counting was performed. The mean value for every detected 
line was determined as the mean of the Gauss distribution. The squares represent 
measured points while the line gives the linear fit to data. Fit parameters are indi- 
cated on the graph. 
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Fig. 13. Fit of measured spectra for oxygen (top), calcium (middle) and tungsten 
(bottom) ions. All spectra are measured with a 40 /is time window. Corresponding 
fit values of the mean number of pe~ per single ion hit are O: 5.15 it 0.07 pe~, 
Ca: 2.78 ± 0.07 pe" and W: 1.81 ± 0.13 pe". 
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Fig. 14. Dependence of the quenching factor of CaW04 on the atomic mass of 
selected elements. The squares are the data of the presented work. The neutron 
scattering result (shown as x) also measured at room temperature is from [16]. The 
circles are from a measurement with a cryogenic detector at a temperature of 7mK. 
The Pb value, (48.7ib7.1), for an energy of ~ 100 keV, is taken from [3], while the 
He value (5.70±0.01) is derivated from the 2.3 MeV a-peak in [3, Fig.7]. 
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